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Abstract 

Using an interface displacement model derived from a microscopic density 
functional theory we investigate thin liquidlike wetting layers adsorbed on flat 
substrates with an embedded chemical heterogeneity forming a stripe. For a 
wide range of effective interface potentials we find first-order phase transitions 
as well as continuous changes between lateral interfacial configurations bound 
to and repelled from the stripe area. We determine phase diagrams and 
discuss the conditions under which these morphological changes arise. 

68.10.-m,82.65.Dp,68.45.Gd 



1 



Typeset using REVTgX 



A variety of experimental techniques have emerged which allow one to endow solid sur- 
faces with a rich, well-defined and permanent chemical pattern while keeping the surface flat 
on a molecular scale (see, e.g., ref. [0). An important application of these structures is mi- 
crofiuidics i.e., the guidance of tiny amounts of adsorbed liquids along these chemical 
microstructures. This enables one to control the microscopic flow of liquids on designated 
chemical channels and to faciliate the fabrication of "chemical chips" which may act as mi- 
crolaboratories for the investigation and processing of rare and valuable liquids |^ . Although 
these applications involve dynamical processes, as a prerequisite it is important to investi- 
gate these systems in thermal equilibrium as a function of the thermodynamic parameters 
pressure (or, equivalently, chemical potential /i) and temperature T. Some recent studies 
were concerned with the static behaviour of liquid channels on chemical lanes within the 
micrometre range (see refs. and 0). On this scale the morphology of the adsorbed liquid 
is determined by gross features such as the various surface tensions involved. However, with 
the rapidly proceeding miniaturisation of microstructures in mind, here we are interested in 
a much smaller scale within which details of the molecular forces become relevant [^. 

The paradigmatic system considered here is a chemically heterogeneous surface which, 
in top view, exhibits a single stripe. The substrate is fiat and composed of two different 
chemical species such that one of them (denoted by forms a single slab of width a 

embedded in the other one (denoted by "— "; see fig. [^(a)). Based on a density functional 
approach p it turns out that the liquid-vapour interface of a thin liquidlike layer in contact 
with the wall interacts with the wall via an effective interface potential A(x, /(x)), where l{x) 
is the local thickness of the layer at the lateral position x. The substrate potential entering 
into A(x, l{x)) has been obtained from a pairwise summation over all substrate- fluid particle 
interactions assuming sharp chemical steps at a: = ±a/2. Since the system we consider here 
is translationally invariant in the y direction, A does not depend on y. We take both the 
substrate-fiuid (s) and the fluid-fluid (/) interaction potential to be of Lennard- Jones type: 
^s,f{f) = 4esj((crs,//r)^^ — (o"sj/r)^). We choose the two chemical species such that a flat, 
semi-infinite, and homogeneous substrate composed of each of the species alone exhibits an 
effective interface potential A±(Z) as depicted in fig. |l](b). (One may interpret these interface 
potentials as corresponding to a "hydrophilic" (A_(/)) and a "hydrophobic" (A_,.(/)) surface 
even for a nonvolatile liquid.) For our choice of potential parameters the outer part of the 
substrate undergoes a critical wetting transition at /c^T^/e/ = 1.2, whereas a homogeneous 
substrate filled with species corresponding to the stripe part exhibits a first-order wetting 
transition at ksT^/ef ~ 1.102. For the temperature ksT/ef = 1.1 considered throughout 
the paper both substrate types are only partially wet at bulk liquid-vapour coexistence 
/i = /io, i.e., Ajj, = jjQ — n = 0. 

Within a simple interface displacement model, which can be derived systematically from 
density functional theory (see ref. |^) the equilibrium contour l{x) of the liquid- vapour 
interface minimizes the functional 



ns[lix)] = j dxdyA{x,l{x))+aigj dx dy ^ 1 + (^^^^^^^ (1) 

with the surface area A = L^Ly of the substrate surface and A(x,/) = AfiApl + uj{x,l) 
where, as it turns out, u{x,l) = X]i>2 '^«(^)^ * ^'^^ Lennard- Jones potentials, Ap = pi — pg 
is the difference in number densities between the bulk phases, and aig is the surface tension 
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associated with the area of the hquid-vapour interface. Instead of the chemical potential one 
may use the pressure p of the bulk vapour phase as a thermodynamic control parameter. 
In this case A/x = corresponds to p = PsatiT) at which the vapour phase is saturated. 
In a pressure-temperature ensemble one has to replace A/xAp by the pressure difference 
Ap = Psat — P- In eq. (|1|) we have omitted the free energy contributions from the wall- 
liquid interface at z = which are constant with respect to l{x). Subtraction of the surface 
free energy VLs{l-) corresponding to the homogeneous outer ("— ") substrate yields the line 
contribution LyVLi[l{x)] = fls[K^)] ~ ^s(^-) to the free energy of the fluid configuration 
associated with the presence of the chemical stripe. f2s[/(a;)] (or equivalently, is 
minimized numerically with respect to l{x) yielding the equilibrium contour within mean 
field theory. 

Using the interface potential A(a;, /) we find equilibrium interfacial morphologies, perti- 
nent examples of which are shown in fig. ^(a). This figure depicts the interface profiles for 
varying stripe widths a at a fixed value A/x = O.OOSej. Within a wide range of values for 
a there are two different minimal interfaces, one "bound" to and the other "repelled" from 
the stripe area. 

Figure |^(b) displays the values of the line free energy density Qi corresponding to the 
interface profiles shown in fig. |^(a). For large a the solution bound to the stripe area has 
a lower line free energy Qi than the solution repelled from the stripe. As a is decreased 
the latter solution is favoured in terms of the free energy. Thus at some value a = at there 
is a phase transition between both interfacial configurations. Due to the break in slope of 
Qi{a) at a = at this transition is first order. We refer to this phenomenon as "morphological 
phase transition" in the sense that the interface profile l{x) undergoes an abrupt structural 
change. 

The repelled solution exhibits a coverage which is even larger than that corresponding to 
the homogeneous "— " substrate. This counterintuitive result shows that compared with a 
homogeneous substrate, one can increase the total adsorption by the immersion of a slab of 
a material that favours thinner liquidlike films. The occurrence of this phenomenon persists 
if the depth of the slab is not macroscopicly large as in fig. |l](a) but only molecularly small 
corresponding to a different material within an imprinted overlayer covering a homogeneous 
underlying substrate. This example shows that gradual changes in the architecture of chem- 
ically microstructured devices can lead to abrupt changes in the morphology of adsorbed 
liquids. We emphasise that such phase transitions are not only of theoretical interest. Their 
existence illustrates the care required to avoid the unwanted filling of the nonwet space 
between liquid channels in such devices. 

For |x| ^ oo the interface profile l{x) asymptotically approaches the equilibrium film 
thickness I- corresponding to the outer part of the substrate. If the stripe width is sufficiently 
large, in the middle of the stripe area l{x) also approaches the equilibrium film thickness 
/+ corresponding to the stripe part. In this expected, l{x) smoothly interpolates 

between the two minima I and II of the interface potentials of the respective homogeneous 
and fiat "+" and "— " substrate taking full advantage of the deep minimum I (see fig. |T](b)). 
If the stripe width shrinks this gain in free energy decreases accordingly and, moreover, the 
relative cost for this benefit in terms of the associated increased area of the liquid-vapour 
interface increases. For a < at the loss of free energy by occupying the higher minimum III 
instead of I is outweighed by the gain in free energy due to a reduced area of the liquid- 
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vapour interface, leading to the repelled solution because the position of the local minimum 
III occurs at a larger value of / than for II. The physical nature of this transition is similar 
to the "unbending" transition on a corrugated wall as described in ref. 0. As explained in 
ref. [|10| similar effects such as "out-of-phase behaviour" can also occur on nonplanar walls 
if there are competing minima of the effective interface potential. 

For values of A/i which are larger than a certain critical value A/ic the morphological 
phase transition does not occur, i.e., there is only one stable solution for every value of a. 
Figure |^ shows the behaviour of the liquid- vapour interfaces for A/i = 0.014e/ > A/ic- In 
this case the interface profile gradually changes from a repelled configuration for small a to a 
bound configuration for large a. The corresponding line free energy Qi is an analytic function 
that does not exhibit a cusp singularity and associated metastable branches. Figure § 
displays the line of first-order phase transitions in the (a, Ayu) plane at which the bound and 
the repelled solutions coexist. Within the present mean field description this line ends at a 
critical point with A/ic ~ O.OlOe/ and Qc ~ 3. 80"/. 

The number of degrees of freedom involved in the morphological phase transition is 
proportional to the volume Ly a (/_ — /+) and thus quasi-onedimensional. Therefore the 
system cannot support a phase transition in the strict sense of statistical mechanics. Actually 
the fluctuations are so strong that at the phase boundary shown in fig. ^ the system cannot 
sustain the sharp coexistence between the bound and the repelled configuration. Instead the 
system will break up into domains along the y direction with alternating regions of increased 
and reduced coverage with the positions of the domain boundaries fluctuating. This leads to 
a rounding of the first-order phase transition due to the aforementioned finite-size effect in 
two directions and eliminates the critical point shown in fig. ^. However, for large values of a 
or /_ — the number of degrees of freedom involved becomes quasi-t wo dimensional so that 
the rounding of the first-order phase transition sharpens up turning it ultimately into a true 
discontinuity. Following the general ideas of finite size scaling of first-order phase transitions 
|TT[| one can estimate the width 25fi of the rounding around the mean field location Afit of 
the phase transition according to the implicit equation 

m,{T, a, AMa) + S^.^a))\ < M exp (^_ -^'fTMy j (2) 

where is the bulk correlation length (which is of the order of o"/) and ^ <^ig{l- — ^+) 
estimates the cost in free energy to maintain an interface between two domains as described 
above; Afi/ is the difference in line free energies between the bound and the repelled solution 
(compare refs. |12| and A rough numerical estimate of 5^ yields the boundaries of the 



smeared out transition region which are indicated by the dashed lines in fig. ^. Whereas far 
from the mean field critical point the width of the transition region is exponentially small 
such that the transition is quasi-first order, for A/x increasing towards A/ic the rounding of 
the phase transition becomes significant. 

In closing we note that the only prerequisite for the occurrence of the morphological 
phase transition described here is that the effective interface potentials of the materials 
involved exhibit the gross features shown in fig. 0(b). This holds both for volatile and 
nonvolatile liquids and does not depend on a particular type of interaction potential or 
thermodynamic ensemble. Such morphological transitions also occur for other shapes of 
chemical substrate heterogeneities such as periodic arrays of stripes and circular or even 
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irregular areas. A variety of experimental techniques such as reflection interference contrast 
microscopy (RICM) |]14| or a force microscope used in tapping mode are being developed 
which allow one to scan the morphology of fluid films with a spatial resolution down to 
fractions of a nanometre. The detection of the morphological phase transition predicted 
here could serve as a promising testing ground for the development of such experimental 
techniques. 
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FIG. 1. (a) Cross section of the morphology of a hquidlike wetting film covering a planar 
substrate that contains a slab of different material. The substrate surface is located at z = and 
the system is translationally invariant in the y direction. Viewed from the top the stripe region is 
filled with chemical species denoted by "+" and extends from x = —a/2 to x = a/2, i.e., the stripe 
width is a. I- = l{\x\ — > oo) is the equilibrium thickness of the liquidlike wetting film corresponding 
to a flat, homogeneous substrate composed of "— " particles alone. The stripe and the outer region 
favour a thin and a thick wetting layer, respectively, (b) The effective interface potentials ^+{1) 
(full line) and A_(Z) (dashed line) of the flat and homogeneous substrates composed of the chemical 
species "+" and "— ", respectively, are minimized for / = and 1 = 1-, respectively. These interface 
potentials correspond to ksT/ef = 1.1 and A///e/ = 0.003 (as, c.f., in fig. §). For the 

positions of the minima II and III are shifted towards larger values of / and to / = oo, respectively, 
whereas I remains basically unchanged. The global minima I and II of A-|-(/) correspond to the 
equilibrium film thickness of the liquidlike film adsorbed on the respective homogeneous substrate. 
A-|-(/ — > oo) increases linearly as AfiApl. A±{Afi = 0,1 = oo) = and A±{Ajj, = 0,1 = hji) < 0, 
i.e., at two-phase coexistence A/i = and at the above temperature both types of substrates are 
only partially wet. 
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FIG. 2. (a) Liquid-vapour interfaces on the stripe-shaped heterogeneity calculated for 
ksT/ef = 1.1 and llfi/ef = 3 • 10^'^ based on the effective interface potentials shown in fig. |l](b). 
The interface profiles are calculated for a/uf = 10 (dashed line), 20 (dotted lines), and 40 (full 
lines). For a = lOo"/ there is only a repelled solution. The dash-dotted lines indicate the equilib- 
rium film thicknesses (— •) and I- (— • •). (b) That part of the line contribution VLi to the grand 
canonical potential which depends functionally on l{x), for the same choice of parameters as in (a). 
At a = at ~ 17 Act f there is a first-order phase transition from the repelled (dashed line) to the 
bound solution (full line). The dot indicates the end of the metastable branch of bound solutions 
for a < at- The metastable branch of repelled solutions extends far to large values of a. 
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FIG. 3. Profiles of the liquid- vapour interface exposed to substrates as shown in fig. |^(a) for 
kBT/ej = 1.1 and A/x = 0.014ej > Afi^- The numbers accompanying each graph indicate the 
corresponding value of the stripe width a in units of cry. The same set of interaction parameters is 
used as for fig. |2|. The dash-dotted lines indicate the equilibrium film thicknesses (— •) and l- 
(— • •) (see fig. |l|(b)). Since the system is above the critical point (A^c,«c) of the morphological 
phase transition there is only one stable solution for every value of a. There is no first-order phase 
transition and 0;(a) is a smooth function. 
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FIG. 4. Line of first-order phase transitions between solutions bound to and repelled from 
the stripe area (full line). The same set of interaction parameters is used as for figs. ^ and ^. 
Within mean field theory the line ends at a critical point which is indicated by the diamond 
symbol. For decreasing Afi the line of phase coexistence extends further up to higher values of a. 
The dashed lines indicate the region within which the rounding of the first-order transition due 
to thermal fluctuations takes place. For decreasing values of a, i.e., upon approaching the mean 
field critical point, this rounding becomes more and more important. Thus the fluctuations erase 
the phase transition near Oc and eliminate the associated critical point. For large values of a the 
fluctuation-induced rounding of the first-order phase transition is practically not detectable. 
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